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Abstract 

Ground-state properties of the Blume-Emery-Griffiths model with antiferromagnetic 
nearest-neighbor interactions on a triangular lattice are investigated in the presence 
of an external magnetic field. In particular, we explore the model's parameter space 
and identify regions with different degenerate ground states that may give rise to 
different magnetic phases also at finite temperatures. We demonstrate the presence 
of such phases by Monte Carlo simulations of magnetization processes for selected 
values of parameters. 
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1 Introduction 



A frustrated triangular lattice Ising antiferromagnet (TLIA) with spin 1/2 is 
long known to display no long-range ordering down to zero temperature [I], albeit 
the ground state is critical with the power-law decaying spin-correlation function [2]. 
Nevertheless, a long-range order can occur in the ground state if the spin is larger 
than some critical value, estimated as 11/2 [3H5]. Generally, the lack of order in 
frustrated spin systems is due to large ground-state degeneracy. However, this can 
be lifted by various perturbations, such as an external magnetic field [6H9] or se- 
lective dilution [IOl|Tl], which can result in long-range ordering even in TLIA with 
spin 1/2. In the Ising models with spin larger than 1/2, a single-ion anisotropy and 
higher-order (e.g., biquadratic) exchange interactions may play a crucial role in their 
critical properties (see, e.g., p^HT^ ). The model that incorporates the above men- 
tioned effects is known as Blume-Emery-Griffiths (BEG) model [H] and has a long 
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history of investigation [Ti]423] . In the case of the BEG model with antiferromagnetic 
interactions, it is interesting to study its behavior in an external magnetic field. A 
recent study by using exact recursion relations on the Bethe lattice ^24j produced 
some interesting results, such as the reentrant phenomenon with the consecutive 
phase transitions from disorder to order and back to disorder as the field was in- 
creased. However, as already observed for example in the simple TLIA model with 
spin 1/2, a frustrated antiferromagnet in the presence of an external field behaves 
quite differently than its nonfrustrated counterpart [6]-|9]. A frustrated Blume-Capel 
model on a triangular lattice, which is special case of our BEG model with zero 
biquadratic exchange interactions and zero magnetic field, has been investigated 
by position-space renormalization group methods [25] and has been shown to dis- 
play finite-temperature long-range ordering within a certain range of the single-ion 
anisotropy strength, accompanied with a multicritical behavior. In order to under- 
stand finite-temperature behavior of the frustrated spin systems, which is often 
unexpected and intricate, it is important to understand their ground-state proper- 
ties. 

In the present Letter we consider the geometrically frustrated antiferromagnetic 
BEG model on a triangular lattice in the presence an external magnetic field. As 
a result of the frustration, the ground state of such a system is highly degenerate. 
On the other hand, the model features a number of different perturbations that can 
lift this degeneracy in different ways and thus we can expect much richer variety of 
magnetic structures in a broad parameter space than for a nonfrustrated case. 



2 Model and methods 

We consider the spin-1 Ising model on a triangular lattice described by the 
Hamiltonian 

n = -JiY. SrS, - J2 E s^s] -DY^sf -h^^s,, (i) 

where = ±1, is a spin on the ith lattice site, denotes the sum over nearest 
neighbors, Ji < is an antiferromagnetic bilinear exchange interaction parameter, 
J2 is a biquadratic exchange interaction parameter, Z) is a single-ion anisotropy 
parameter, and h is an external magnetic field. 

Considering the presence of only nearest-neighbor interactions, it is sufficient to 
focus on an elementary triangular plaquette formed by the neighboring spins. Then, 
a reduced zero-temperature energy per spin can be expressed as 

e/l Jil = E S,S, - J2/I Jil E SlSf - D/{3\J,\)Y.Sl - V(3| Ji|) E^^> (2) 

{k,l) {k,l) k k 

where the summation {k, I) runs over the nearest neighbors Sk and Si {k, I = 1,2, 3) 
on the plaquette. Ground-state (GS) spin configurations Sgs = {5"!, S'2, S'3} can be 
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found as configuration^^ that minimize the energy functional /(S, P) = e/| Ji| with 
the parameters values P = (J2/I Ji|, -D/| Ji|, Ji|), i.e.: 

SGs = argmin/(S,P). (3) 
s 

Apparently, for some parameters values there will be no single optimal solution but 
rather a multitude of degenerate states with the same energy. Then the GS energy 
per spin can be obtained as a function of the parameters P: /gs(P) = /(Sgs? P)- 

In order to demonstrate the presence and character of different phases resulting 
from the GS configurations found from Eq. ([3]), we further employ Monte Carlo 
(MC) method and study magnetization processes in different regions of the param- 
eter space. We perform MC simulations on a spin system of a moderate linear size 
of L = 24, employing the Metropolis dynamics and applying the periodic bound- 
ary conditions. For thermal averaging we consider = 5 x 10^ MCS (Monte Carlo 
sweeps or steps per spin) after discarding another 10^ MCS for thermalization. For 
configurational averaging and calculating error bars we carry out five independent 
runs. The simulations are performed at a fixed reduced temperature close to the 
ground state {t = kBT/\Ji\ = 0.1) and a varying external field h/\Ji\. The simula- 
tions start from zero field, using a random initial configuration, and then the field 
is gradually increased with the step Ah/\Ji\ =0.1 and the simulation starts from 
the final configuration obtained at the previous field value. The triangular lattice is 
considered to be consisting of three interpenetrating sublattices A, B and C, and we 
calculate the respective sublattice magnetizations per site 




and the total magnetization per site 

^ = ( E Sn)/L', (5) 

^ n=l ' 

where (■ ■ ■ ) denotes thermal averages. 



3 Results 

By minimizing the energy functional (Eq. [3]) we identified GS configurations and 
calculated GS energy in a broad parameter space. In Fig. [1] we show the GS energy 
surfaces in the D/\Ji \ — h/\Ji \ parameter plane for various values of the biquadratic 
exchange parameter J2/\Ji\- Superimposed are the lines separating different phases, 
which are presented in Table [H with the lists of the corresponding degenerate states 

^ There are in total 27 configurations to be considered. 
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and their energie^. The analytical expressions for these phase boundaries can be 
obtained by pairwise equating the GS energies. We can see that, depending on the 
parameters values, the system can be in up to eighth different phases. In the asymp- 
totic limits of the parameters, for large negative values of the single-ion anisotropy 
Dl\Ji\ the system is in the nonmagnetic state Sqs = {0,0,0}, while for large posi- 
tive values of D /\Ji\ the nonmagnetic spin states are suppressed and we can expect 
the behavior similar to that of the spin-1/2 model. On the other hand, for large 
magnetic fields h/\Ji\ the system adopts the magnetic state Sgs = {1; 1; 1} with all 
the spins fully aligned with the field direction. In the mentioned cases of the large 
negative D/\Ji\ ((T)) and large h/\Ji\ ((t)) the degeneracy is lifted completely. On 
the other hand, there is still three-fold ((2),(3),(6}) and six-fold ((4),(5),(8}) degen- 
eracy if the parameters are of intermediate and relatively small values, respectively. 
We note that at the boundaries the degeneracy is further increased due to contribu- 
tions of states from the neighboring phases and there are up to five points in which 
even three different phases coincide. 

[Fig. 1 about here.] 

[Table 1 about here.] 

With the respect to the above arguments on the behavior of the local spin pat- 
terns, we can expect the presence of interesting magnetic structures on the entire 
lattice in different regions of the parameter space. The character of the respective 
phases can be revealed by inspecting the sublattice magnetizations mA, m-Q and mc 
as well as the total magnetization m, obtained from MC simulations. Their field- 
dependencies for a fixed temperature t = 0.1 and the biquadratic exchange param- 
eter J2/\Ji\ = 0, and three different values of the single-ion anisotropy parameter 
are plotted in Fig. |2J The inset spin snapshots in the respective figures illustrate 
typical spin configurations for different values of the field h/\Ji\, representing dif- 
ferent phases. Positive values of the parameter D tend to suppress nonmagnetic 
states and the spin system is expected to behave like a spin-1/2 TLIA. Namely, it 
is expected to show no long-range order in zero field, the ferrimagnetic state of the 
type {mA,mB,mc) = (1, —1, 1) with two sublattice magnetizations aligned parallel 
and one antiparallel to the field direction (i.e., 1/3 magnetization plateau) within 
h/\Ji\ G (0,6), and the ferromagnetic state with all three sublattice magnetizations 
aligned with the field direction for h/\Ji\ > 6 [6H9]. This scenario is supported by 



the plots presented in Fig. 2(a) for D/\Ji\ = 1. On the other hand, negative values 



of the parameter D tend to enhance nonmagnetic states. However, for relatively 
small magnitudes, more specifically for J2/IJ1I = it is within —3/2 < D/\Ji\ < 



hence also for D/\Ji\ = —1 shown in Fig. 2(b), in zero field the minimum energy 
is realized if nonmagnetic states are restricted to only one sublattice, while the 
remaining two sublattices are ordered antiferromagnetically, i.e., the structure of 
the type (1,-1,0). As the field increases, the structure evolves first to the type 



The phases (4) and (s) are present only at h/\Ji\ = 0. 
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(1, —1, 1), then to (1, 0, 1) and finally for sufficiently large fields to (1, 1, 1). Thus, 
the total magnetization m features four plateaus with the values mi = 0, m2 = 1/3, 
77^3 = 2/3 and = 1, corresponding to the respective magnetic structures. Finally, 
for negative D of large magnitudes, namely for J2/\Ji\ = it is for D/\Ji\ < —3/2 



and thus also for D/\Ji\ = —2 shown in Fig. 2(c), the nonmagnetic state (0,0,0) 
is the ground state at small enough fields. The increasing field then brings the sys- 
tem successively to the states of the types (1, 0, 0), (1, 0, 1) and (1, 1, 1). This again 
splits the magnetization curve m into four plateaus of the same heights as for the 
case D/\Ji\ = —1, albeit now the steps correspond to generally different magnetic 
structures. 

In Fig. |3]we present similar plots also for a non-zero value of the biquadratic ex- 
change interaction parameter, namely for J2/\Ji\ = —0.5, and three different values 
of the single-ion anisotropy parameter D/\J\ = —1,0, and 1. The phases involved 
are the same as for the case of J2/\Ji\ = and, therefore, we do not illustrate their 
nature again by showing the snapshots. Nevertheless, the sequence of the phases 
through which the system passes as the field increases can be different, such as 
(5) — 7- (2) —7- (3) — !■ (7) for D/\J\ = 0, which cannot be observed at any value of 
D/\J\ if Js/I^il = 0. 

[Fig. 2 about here.] 
[Fig. 3 about here.] 



4 Conclusions 



We have shown that the geometrically frustrated BEG model in an external 
magnetic field displays a relatively rich variety of different ground-state phases in 
the parameter space. They result from the large degeneracy, which is lifted in several 
stages with the increasing influence of the respective parameters. For selected pa- 
rameter values we confirmed the existence of multiple magnetic structures, reflected 
in the step-wise field-dependence of the magnetization curves, in the Monte Carlo 
simulations at very low temperature. We believe that the present results will stimu- 
late further study of the effects of thermal excitations on the individual ground-state 
spin arrangements and to establish entire finite-temperature phase diagrams. 
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List of Figures 



1 (Color online) Surfaces of the ground-state energy /qs in the 

D /\Ji\ — h/\Ji \ parameter plane for various values of the biquadratic 
exchange parameter J2/\Ji\- The lines mark the borders between 
different phases, presented in Table [H 



(Color online) Field-dependencies of the total magnetization m 

and sublattice magnetizations mA, m-Q and mc at the temperature 

t = 0.1 for J2/\Ji\ = and (a) D/\Ji\ = 1, (b) D/\Ji\ = -1 and 

(c) D/\Ji\ = —2. The snapshots in the respective figures illustrate 

typical spin configurations in different phases. The red (light gray), 

blue (dark gray) and white circles represent the spin states —1, +1 

and 0, respectively. 

(Color online) Field-dependencies of the total magnetization m 

and sublattice magnetizations m\, me and mc at the temperature 

t = 0.1 for J2/IJ1I = -0.5 and (a) D/\Ji\ = 1, (b) D/\Ji\ = and 

(c) D/lJil = -1. [10 
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Fig. 2. (Color online) Field-dependencies of the total magnetization m and sublattice 
magnetizations mA, w-b and mc at the temperature t = 0.1 for J2/\Ji\ = and (a) 
D/\Ji\ = 1, (b) D/\Ji\ = —1 and (c) D/\Ji\ = —2. The snapshots in the respective figures 
illustrate typical spin configurations in different phases. The red (light gray), blue (dark 
gray) and white circles represent the spin states —1, +1 and 0, respectively. 
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Fig. 3. (Color online) Field-dependencies of the total magnetization m and sublattice 
magnetizations mA, w-b and mc at the temperature t = 0.1 for J2/\Ji\ = —0.5 and (a) 
D/\Ji\ = 1, (b) D/\Ji\ = and (c) D/\Ji\ = -1. 
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List of Tables 



Ground-state phases identified by Eq. ([3]) witli tlie corresponding 
degenerate states Sgs and energies /gs(P), at the parameters values 

P=(J2/|Jl|,D/|Jl|,VI'/l|)- 
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Table 1 

Ground-state phases identified by Eq. ^ with the corresponding degenerate states Sgs 
and energies /gs(P)) at the parameters values P = {J2/\Ji\, D/\Ji\, h/\Ji\). 

Sgs = {'S'l, 5*2, S3} 



/gs(P) 



{0,0,0} 
{0,0,1}; {0,1,0}; {1,0,0} 
{0,1,1}; {1,0,1}; {1,1,0} 
{-1,0,0}; {0,-1,0}; {0,0,-1}; {0,0,1}; {0,1,0}; {1,0,0} 
{-1,0,1}; {-1,1,0}; {0,-1,1}; {0,1,-1}; {1,-1,0}; {1,0,-1} 
{-1,1,1}; {1,-1,1}; {1,1,-1} 
{1,1,1} 

{-1,-1,1}; {-1,1,-1}; {-1,1,1}; {1,-1,-1}; {1,-1,1}; {1,1,-1} 
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